25 September 1997, to be published in PHYSICA E 



Semiclassical Approach to Orbital Magnetism of 
Interacting Diffusive Quantum Systems 

D. Ullmo, (1 ' 2) K. Richter/ 3 ) H.U. Baranger/ 1 ) F. von Oppen/ 4 ) and R.A. Jalabert (5) 

W Bell Laboratories-Lucent Technologies, 700 Mountain Ave., Murray Hill NJ 07974 
^Division de Physique Theorique, Institut de Physique Nucleaire, 91406 Orsay Cedex, France 
( 3 ) Max- Planck- Institut fur Physik komplexer Systeme, Nothnitzer Str. 38, 01187 Dresden Germany 

^Department of Condensed Matter Physics, Weizmann Institute, 76100 Rehovot, Israel 
(5) Universite Louis Pasteur, IPCMS-GEMME, 23 rue du Loess, 67037 Strasbourg Cedex, France 

Abstract 

We study interaction effects on the orbital magnetism of diffusive mesoscopic 
quantum systems. By combining many-body perturbation theory with semi- 
classical techniques, we show that the interaction contribution to the ensemble 
averaged quantum thermodynamic potential can be reduced to an essentially 
classical operator. We compute the magnetic response of disordered rings and 
dots for diffusive classical dynamics. Our semiclassical approach reproduces 
the results of previous diagrammatic quantum calculations. 



A. Introduction 

The interplay of disorder and interactions in mesoscopic systems has attracted consider- 
able attention Interaction effects on transport through small quantum dots || as well as 
on thermodynamic properties like persistent currents and orbital magnetism are of present 
interest. In the latter case, the unexpectedly large measured persistent current of small 
metal rings [[I]-[5|] pointed towards the importance of such interaction effects and motivated 
a large number of theoretical approaches ||. 

For the description of thermodynamic quantities, semiclassical expansions have proven 
particularly useful, both within the independent-particle model and for interaction 

effects p3l , |TTll . These studies established a close relation between the classical dynamics 
and the quantum-mechanical magnetic response. In particular, studies of ballistic systems 
showed that the quantum thermodynamic properties are sensitive to whether the classical 
dynamics is regular or chaotic [p|,p|,|TT| . 

In this paper we apply these semiclassical techniques to the orbital magnetism of inter- 
acting systems whose non-interacting classical dynamics is diffusive. Specifically, we present 
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semiclassical derivations of the interaction contributions to the persistent current of metal 
rings and to the susceptibility of singly-connected two-dimensional diffusive systems. We 
recover results obtained previously by quantum diagrammatic calculations pT2|-|Hj1, showing 
that the semiclassical approach is on the same level of approximation. By semiclassically 
evaluating the relevant diagrams appearing in the many-body perturbation series for the 
thermodynamic potential, we express the latter in terms of an essentially classical operator. 
This expression provides a convenient starting point for further calculations. Moreover, by 
making the connection with the classical dynamics, it provides a physically intuitive picture 
of the interplay between disorder and interaction. 



B. Diagrammatic perturbation theory 

We are interested in the orbital magnetism of a mesoscopic quantum system subject to 
an external magnetic field B. While the magnetic response of a singly-connected system 
is usually measured in terms of its susceptibility \i t ne magnetic moment of a ring-type 
structure threaded by a flux = BA (where A is the enclosed area) is usually described by 
the related persistent current /. Both are given in terms of the thermodynamic potential Q 
as (V being the area (volume) of the structure) 

' X = ~VdB^ () 
To calculate the interaction contribution to the magnetic response, the high-density 



expansion (RPA) of the thermodynamic potential [pl| has to be extended by including in- 
teraction corrections from diagrams with the Cooper-channel. This was originally performed 
in the context of superconducting fluctuations and then applied to disordered normal metals 
P^T7[. Such expansions usually yield reliable results even beyond the high density limit, 
if the relevant sets of terms are properly resummed. The relevant Cooper-like diagrams are 
shown in Fig. |l|. The screened Coulomb interaction (wavy lines) can be treated as local 
EUH: U(r - r') = A A^(0)" 1 5(r - r'). Here, JV(0) denotes the density of states and the 
bookkeeping index Ao = 1 identifies the order of perturbation. For the local interaction, 
direct and exchange term are equivalent up to a factor of (—2) due to the spin sums and the 
different number of fermion loops. The corresponding perturbation expansion for this inter- 
action contribution Q to the thermodynamic potential, which yields the magnetic response, 
can be formally expressed as 
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FIG. 1. Leading Cooper-channel diagrams for the interaction contribution to the thermody- 
namic potential. 
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Here, a; denotes the bosonic Matsubara frequencies uo = 2nir/f] with (3 = 1/kT. The particle- 
particle propagator £(u;) is expressed (in position representation) in terms of products of 
finite-temperature Green functions as |T8| 

S r,r' (W) = ^yTg) E ^r,r' (e)C?r,r' (w - g) . (3) 

Here, the sum runs over the fermionic Matsubara frequencies e = (2n + 1)tc//3. The short- 
length (high-frequency) behavior is included in the screened interaction, thus requiring a 
cutoff of the frequency sums at the Fermi energy E F fl4j| . The straight lines in Fig. 1 
represent finite-temperature Green functions of the non-interacting system. They are of 
the form 

g r>T ,{e) = 9(e)G* r ,(E F + ie) + 6(-e)G* r ,(E F +ie) (4) 

in terms of the retarded and advanced Green functions G R,A which are related by G^ r , (E) = 
[G$ T (E*)]*. For diffusive systems, they include the presence of the disorder potential. 



C. Semiclassical formalism 

Both in ballistic and diffusive samples, the Fermi wavelength \p is often the shortest 
lengthscale. It is in this situation that we can apply semiclassical techniques to compute 
H rr t(u). Here, we will moreover assume that the magnetic field B is classically weak, i.e., 
that the cyclotron radius R c 3> min{Z, L} (with I the elastic mean free path and L the 
system size). 

Semiclassically, the retarded Green function is represented as a sum of contributions G^'J, 
over all classical paths j from r to r' ||19|| , 

G?A E )~ E Dje^-^ 2 . (5) 

j-.r—ir' 

Here Sj— J r r p • dr is the classical action of trajectory j. The prefactor Dj includes the 

classical phase space density [Dj = (1/ ^27r(ih) 3 xx') \d 2 Sj/ dydy'^ 2 in two dimensions]. Vj is 
a Maslov index. The semiclassical approximation makes the temperature and magnetic-field 
dependences of the finite-temperature Green function transparent. Employing (dSj/dE)=tj 
and (dSj/dB)=(e/c)Aj, where tj and Aj are the traversal time and area, one finds 

G?$(E F +ie, B) ~ G^,{E F ,B = 0) x exp [-etj/H] x exp foirBAj/fa] (6) 

where 0o — he/ e is the flux quantum. Note that temperature exponentially suppresses the 
contributions of long paths to each Green function. 
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Semiclassically, the particle-particle propagator Ti T)T i(u) is then represented as a sum over 
pairs of paths between r and r'. Off-diagonal pairs (of different paths) generally contain 
highly oscillatory contributions which do not survive an ensemble (disorder) average. (There 
can be exceptions as discussed in Ref. [IT||.) On the other hand, the diagonal pairing of 



each orbit j with its time reverse persists upon averaging since their dynamical phases 
exp[iSj(B = 0)/h] cancel while retaining a magnetic-field dependence. A more detailed 
semiclassical analysis [TT|] shows that the Cooper series in Fig. 1 contains the magnetic-field 
sensitive contribution to Q which is leading order in h. 

Using Eqs. (f|), (|5|) and (^) in Eq. (^) and performing the Matsubara sum yields for the 
diagonal part of £ 
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(7) 



The sum runs over all trajectories longer than the cutoff Ao = A F /7r [corresponding to the 
upper bound E F on the Matsubara sum in Eq. The temperature dependence in Eq. (^) 
enters through the function R(x) = xj sinh(x) introducing the time scale 



hp 

t T = — 
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(8) 



and the related length scale Lt = v F tT, with v F being the Fermi velocity. This semiclassical 
framework allows us to reduce the original quantum problem to E^ which no longer exhibits 
variations on the quantum scale but only on classical scales. We emphasize that the 
representation (0) of E^ is rather general since we have not yet made any assumption 
about the classical dynamics of the system. In particular, it applies to both diffusive and 
ballistic systems. On the basis of Eq. (|7]), we have recently studied interaction effects 
in ballistic quantum dots 



11 



Specifically, we show that the interaction induced orbital 
magnetism scales differently for systems with regular and chaotic non-interacting classical 
counterparts. 

Here, we focus on diffusive systems for which it is useful to relate YS D > to classical 
probabilities satisfying the diffusion equation. To this end we introduce an additional time 
integration in Eq. (0) and make use of the relation [^] 



2tt 2 



N(0) 
2nh 



P(r,r';t) 



(9) 



between the weights \Dj\ 2 and the classical probability P(r,r';t) to propagate from r to r' 
in time t. 

An n-th order contribution to Q in Eq. (fj) then contains expressions for the joint return 
probability P(ri, . . . , r n , ri; t%, . . . , t n \A) to visit the points (with U being the time between 
Yi and rj + i) under the condition that the enclosed area is A. In diffusive systems, the proba- 
bility is multiplicative, namely / dri . . . dr n P(ri, . . . , r n , ri; ti, . . . , t n \A) = J drP(r, r; ttot\A) 
with ttot = The contribution to Q in Eq. (pj) from the diagonal terms E*-™ then yields 
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n(D) = T, Q i D) = ^J dr J dtcoth (J-) K(t)A(r,t;B) 



(10) 



where the coth(t/tr) arises from the o;-sum in Eq. (^]) and 
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A(r, t)B) = J dA cos ( j P(r, r; t\A) . 



(11) 
(12) 



K(t) accounts for temperature effects while A contains the field dependence and the classical 
return probability. Eqs. (|i"0|)-(|12|) are a general and convenient starting point to compute 
the orbital response of disordered systems. 



D. Diffusive rings 

We start with the computation of the first-order interaction contribution, fl[ , to il- 
lustrate the main ideas. Consider a (thin) disordered ring of width b, crossection a and 
circumference L. For L 3> l,b the motion of particles around the ring effectively follows a 
law for one-dimensional diffusion. Since the area enclosed is given in terms of the number 
m of windings around the ring, one has 

^11 / rn 2 L 2 \ p / t mL 2 \ , . 



m=— oo 



where -D = v F l/d is the diffusion constant (in d dimensions). Because of the disorder average 
the classical return probability does not depend on r. In first order we have 

K 1 (t)=\ R(2t/t T )/t. (14) 

Combining this with the coth function in Eq. ([It]) we find 

^1 = A 2^ COS 7 0m CO (15) 

^ m=-oo \ 90 J 

with 

^ r ^ gCW ex P [-(mL) 2 /(4Dt)] 
9m[1) ~Jo t 2 ^Di • { ' 

Taking the derivative with respect to the flux, we recover the first-order interaction contri- 
bution to the persistent current, first obtained in [15] by purely diagrammatic techniques, 

/i = A 2^ msm — 7 — \9m\T). (17) 

T m=-oo \ <P0 J 

Semiclassically, this first order result was already derived by Montambaux |10| . 
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In addition, our semiclassical approach allows us to obtain the renormalization of the 



coupling constant |12Hl6| due to the higher-order diagrams of the Cooper series. Including 
these diagrams amounts to using the full kernel K(t) in Eq. ([It]) instead of K\{t). Introducing 
the Laplace transform of Ki(t), 

n F -I 

/(PH4A °£ ^ + 2(2, i + 1 ) <18) 
(jif=(3Ef/27t= k F Lr/4), K(t) is given by the inverse Laplace transform 

1 r+ioo A 1 

K(t) = — / dpe** ln[l + f(p)] ~ ————K^t) ; L* = wm(y F t, L T /A) . (19) 
2m J-ioo A m{k F L*) 

The last equality is valid when In k F L* ^> 1 which is certainly satisfied when In k F l ^> 1. 
Therefore, the higher-order terms merely lead to a renormalization of the coupling constant, 
thus reducing the predicted magnitude of the persistent current. In the high temperature 
regime (Lt <C L m ) the coupling constant is renormalized to l/ln(fc F Lr/4). Introducing 
L m =v F (mL) 2 /AD, the average length of a trajectory diffusing m times around the ring, one 
gets at low temperature (L T 3> L m ) a replacement of A = 1 by 1/ \n(k F L m ). These two 
limits agree with results obtained diagrammatically by Eckern Jl6| . 



We note that the functional form of the temperature dependence (exponential T— 
damping |15[) is in line with experiments 0-0] while the amplitude of the persistent current 
with renormalized coupling constant is smaller than the experiments by a factor of ~ 5. 



E. Diffusive two dimensional systems 

Contrary to rings, the geometry imposes no shortest length for returning paths in singly- 
connected systems. One therefore expects a different temperature dependence of the mag- 
netic response. 

Consider a two-dimensional singly-connected diffusive quantum dot. In view of the 
general renormalization property of diffusive systems Eq. (|T9"|), the diagonal part of the 
thermodynamic potential from the entire Cooper series [Eq. flip]) ] can be written as 

fiW = - fdrfdt — % R 2 (-) A(r,t;B) . (20) 



/3 J J \n(k F v F t) t 2 \t T 

Here we have used L* = v F t in flT9"| ) since the R 2 factor ensures that the main contribution 
to the integral comes from t < tx- In two dimensions the conditional return probability, 
entering into A, is conveniently expressed in terms of the Fourier transform |7J 

P(T,r,t\A) = ± [dk\k\e ikA e M-\k\Dt) 
K ' ' 1 ; 4vr 2 J 11 1 - exp(-2|Jfe|Dt) V ; 

from which one obtains 
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Here we introduced the magnetic time 



J 2 



t B = — — = . (23) 

It is related to the square of the magnetic length L\ which denotes the area enclosing one 
flux quantum (assuming diffusive dynamics). Note that the function R in Eq. (p2|) has a 
different origin than in Eq. (J20|). 

Using the expression fl22"|) in Eq. (|20|) and taking the second derivative with respect to 



the field, we find for the susceptibility, 

X (D) = _ 6 [°° dt 2 (t_\ (t_\ 

\Xh\ X F ' Jr el t\ll{k F V F t) \t T ) \t B ) K ' 

where R" is the second derivative of R. The susceptibility is normalized to the two- 
dimensional diamagnetic Landau susceptibility xl = —e 2 /(12irmc 2 ). 

In the above time integral the elastic scattering time T e \ = l/v F enters as a lower bound. 
This cutoff must be introduced since for backscattered paths with times shorter than r e \ 
the diffusion approximation ( PT| ) no longer holds [|20|| . On the other hand Eq. ( ]2"4T ) holds 
true only as long as the upper cutoff time t* = min^r,^) is smaller than the Thouless 
time t c = L 2 /D (with L being the system size). For times larger than t c the dynamics 
begins to behave ergodically, and the two-dimensional diffusion approximation is no longer 
valid. Assuming t* < t c , Eq. (|24l) can be approximately evaluated by replacing R{t/tT) and 
R"(t/tB) by -R(O) = 1 and R"(0) = —1/3, respectively, and introducing the upper cutoff t* 
in the integral. The remaining integral yields for t* ^> r e i 

dt - In ( ln ^ F Vp min ( tr ' 1 (25) 



Bl t\n(k F v F t) { \n(k F l) 

The log- log form produced by the 1/tlnt dependence results from the wide distribution 
of path-lengths in the system- there are flux-enclosing paths with lengths ranging from 
about v F r e i up to v F t*. In contrast, in the ring geometry discussed in the previous section 
the temperature dependence is exponential because the minimum length of flux-enclosing 
trajectories is the circumference. 

The averaged susceptibility of a diffusive two-dimensional structure then reads 

IxlI tt { ln(M) J 

One thus finds a log- log temperature dependence for t-r < t B and a log- log B dependence 
for tx > t B . With regard to magnitude, the magnetic response of diffusive systems is 
paramagnetic and enhanced by a factor k F l compared to the clean Landau susceptibility xl- 
Eq. ( f2~6|) agrees with results from Aslamazov and Larkin fll2" |, Altshuler, Aronov and 
Zyuzin [rjJTJL and Oh, Zyuzin and Serota |HJ obtained with quantum diagrammatic per- 
turbation theory. The equivalence between the semiclassical and quantum approaches to 
diffusive systems may be traced back to the fact that the "quantum" diagrammatic per- 
turbation theory relies on the use of the small parameter l/k F l which can be viewed as a 
semiclassical approximation. 
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F. Conclusion 



To conclude, we developed a semiclassical approach to evaluate the interaction contribu- 
tion of the grand potential in a high density perturbative expansion. We showed that the 
averaged quantum magnetic response can be expressed in terms of an operator containing 
the classical probability for particles to return. As an application we computed the orbital 
magnetic response of diffusive rings and two-dimensional quantum dots arising from the 
combined effects of disorder and interaction. 

RAJ and KR acknowledge support from the French-German program PROCOPE. The 
Division de Physique Theorique is "Unite de recherche des Universites Paris 11 et Paris 6 
associee au C.N.R.S.". 
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